function logitprobMAT = logitprob(VMAT, noise, wbar) 
% Computes logit probabilities based on payoff and noise parameter

onesvec = ones(size(VMAT,1),1);
nomiMAT  = exp( (VMAT - onesvec*max(VMAT))./wbar /noise );  %./ wbar
denomMAT = onesvec*sum(nomiMAT);

logitprobMAT = nomiMAT./denomMAT;